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ABSTRACT 


Measurements  of  cross-structure  impact  pressures  and 
relative  motions  between  the  free  surface  and  the  structure 
at  the  point  of  impact  on  a  1/22  scale  SWATH  T-AGOS  model 
were  analyzed  to  determine  whether  these  parameters  could 
lead  to  procedures  for  estimating  the  pressures  from  ship 
motion  calculations.  Although,  the  impact  pressure,  in 
general,  increased  with  the  square  of  the  relative  velocity 
the  data  showed  considerable  scatter  when  plotted  in  this 
format.  Including  the  effects  of  the  angle  of  the  free 
surface  at  the  point  of  impact  in  the  analysis  did  not 
improve  the  correlation.  If  it  is  assumed  that  the  rms 
pressure  is  proportional  to  the  mean  squared  relative 
velocity  and  that  the  probability  distribution  is 
exponential,  then  the  slamming  characteristic  of  the  SWATH 
ship  can  be  estimated  from  ship  motion  computations.  The 
former  assumption  requires  further  experimental 
verification .  - 

ADMINISTRATIVE  INFORMATION 

This  work  was  performed  under  the  Naval  Sea  Systems  Command  General 
Hydrodynamic  Research  Program  administered  by  the  DTNSRDC  Ship  Performance 
Department.  Funding  was  provided  under  Program  Element  61153N,  Task  Area 
SR0230101,  and  Work  Unit  1562-500  in  FY  86. 

INTRODUCTION 

The  slamming  loads  acting  on  the  cross-structure  of  a  SWA1H  ship  can 
present  a  serious  problem  for  operations  in  a  heavy  seaway.  It  is  therefore, 
important  for  the  designer  to  have  available  the  means  for  assessing  the 
slamming  characteristics  of  a  SWATH  ship  early  in  the  design  stage.  The 
methods  employed  for  monohulls  do  not  appear  to  apply  directly  to  the  SWATH 
ship  (at  least  without  some  modifications)  ,  and  an  investigation  was  made  to 


determine  procedures  for  estimating  slamming  loads  more  appropriate  for  the 
SWATH . 

The  method  developed  for  monohulls  by  Ochi1  assumes  that  a  3lam  occurs 
when  the  point  of  impact  contacts  the  water  with  a  velocity  that  exceeds  a 
certain  critical  value.  An  analytical  expression  can  be  derived  defining  this 
event  by  assuming  that  the  relative  motion  between  the  point  of  contact  and 
the  free  surface  is  a  stationary  random  Gaussian  process.  The  critical  or 
thresh--. id  velocity  for  a  520  foot  (158  meters)  Mariner  as  defined  by  Ocni  from 
model  experiments  is  12  ft/sec  (3.7  m/sec).  Thi3  value  is  Froude  scaled  to 
obtain  the  equivalent  values  for  ships  of  different  lengths.  Ochi  also 
provides  analytical  expressions  for  estimating  extreme  values  of  im.pact 
pressure  based  upon  the  assumption  that  the  relative  velocity  is  a  stationary 
random  Gaussian  process  and  that  the  impact  pressure  is  proportional  to 
velocity  squared. 

The  applicability  of  these  relationships,  developed  for  the  keel  slamming 
of  a  monohull,  to  the  cross-structure  slamming  of  a  SWATH  ship  has  not  been 
verified.  Experimental  data  from  three  dimensional  drop  test  of  wedges  by 
Chuang  and  subsequently  applied  to  a  catamaran  cross-structure  suggest  that 
the  slam  pressure  of  a  SWATH  cross-structure  is  dependent  upon  both  the 
relative  velocity  and  the  impact  angle  at  the  point  of  impact.  Experiments 
were  necessary  to  establish  whether  a  relationship  existed  between  slam 
pressure  and  the  relative  velocity  and  angle  of  impact  between  the 
cross-structure  and  free  surface  upon  impact  for  the  SWATH  ship.  Fortunately, 
a  comprehensive  model  test  program  had  been  initiated  for  support  of  the 
T-AGOS  design  which  included  impact  pressure  measurements  on  the 
cross-structure.  Additional  instrumentation  was  added  at  one  gage  location  to 
measure  cross-structure  relative  velocity  and  angle  of  impact  with  respect  to 
the  free  surface. 


DESCRIPTION  OF  EXPERIMENT  AND  INSTRUMENTATION 

As  previously  indicated,  data  for  these  studies  were  obtained  during  a 
much  broader  model  experiment  program  for  the  T-AGOS  SWATH  design.  A  1/22 
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scale  model  was  instrumented  to  obtain  a  comprehensive  seakeeping  evaluation 
including  slam  pressures  on  the  cross-structure.  All  that  was  needed  for 
these  studies  was  the  addition  of  suitable  instrumentation  for  measuring 
velocity  and  angle  relative  to  the  free  surface  at  the  point  of  impact. 
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Figure  1  -  Location  of  Pressure  Gage 

The  pressure  gage  consisted  of  a  strain-gaged  panel  which  was  calibrated 
to  measure  the  average  pressure  over  an  area  corresponding  to  4  ft.  by  4  ft. 
(full  3cale) .  The  location  of  this  panel  on  the  cross-structure  is  shown  in 
Figure  1.  Three  wire  probes  for  measuring  relative  motion  between  the  free 
surface  and  the  structure  were  mounted  in  a  triangular  pattern  encompassing 
the  pressure  gage.  Figure  2  is  a  sketch  showing  the  probe  spacing  around  the 
pressure  panel.  The  relative  motion  data  were  filtered  to  prevent  aliasing 
and  recorded  on  digital  tape  at  a  rate  of  6  samples  per  second.  This  is  the 
normal  procedure  for  recording  relatively  low  frequency  seakeeping  data. 
Because  of  the  rapid  rise  times  associated  with  a  slam,  which  would  require  an 
extremely  high  sample  rate  to  resolve,  the  pressure  data  were  recorded  on 
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analog  tape  along  with  other  structural  load  data.  A  relative  motion 
measurement  was  included  on  the  analog  tape  to  time  correlate  slam  pressure 
with  the  relative  velocity  and  angle  on  the  digital  tape. 


Figure  2  -  Sketch  of  Relative  Motion  Probes 


The  primary  experimental  program  investigated  the  seakeeping  performance 
of  the  T-AGOS  in  various  sea  conditions  at  several  speeds  and  headings; 
however,  only  two  runs  were  selected  for  this  investigation  based  upon  sample 
size  of  the  number  of  impacts.  Operation  in  head  seas  at  8  knots  in  a  sea 
state  9  and  at  zero  speed  in  a  sea  state  corresponding  to  hurricane  Camille 
produced  a  significant  number  of  impacts  on  the  cross-structure  which  made 
these  data  amenable  to  a  comprehensive  statistical  analysis.  A  much  lesser 
number  of  impacts  occurred  at  other  conditions  which  may  be  relevant  to  the 
T-AGOS  seakeeping  assessment,  but  were  not  considered  to  have  sufficient 
sample  size  to  warrant  detailed  statistical  examination. 


METHOD  OF  ANALYSIS 


Figure  3  is  a  graphic  representation  of  the  digital  data  recorded  by  the 


relative  wave  height  probe  as  a  function  of  time. 
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Figure  3  -  Relative  Motion  Time  History 


The  data  show  the  saturation  of  the  probe  when  deck  wetness  occurs  which  is 
typical  during  a  slam  event.  The  synthesized  time  history  of  the  relative 
wave  height  from  the  digital  data  was  readily  correlated  with  the 
corresponding  channel  on  the  analog  tape  containing  the  impact  pressure  data. 
Thus  each  cross-structure  impact  recorded  on  the  analog  tape  could  be  closely 
associated  with  the  corresponding  digital  time  histories  of  the  relative  wave 
height  probes. 


Relative  velocity  was  computed  by  numerically  differentiating  the 
relative  wave  height  displacement.  Three  points  were  used  in  a  polynomial  fit 
to  find  the  value  of  the  derivative  at  the  end  point  which  is  given  by 
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a t(t,)  -d.  [£<«„)  -  +  3f(t2>l 

where,  h  -  time  interval  between  data  points. 

The  values  of  the  computed  velocity  from  the  three  relative  wave  height  probes 
were  averaged  to  obtain  the  velocity  at  the  center  of  the  pressure  gage.  The 
impact  was  assumed  to  occur  at  the  instance  prior  to  deck  wetness  as  indicated 


by  the  flattening  cut  of  the  relative  wave  heig.ot  time  history.  (See  figure 
3  .  ) 

The  angle  of  impact  was  calculated  from  the  relative  wave  height 
measurements  assuming  the  free  surface  to  be  a  simple  plane  surface  between 
the  three  probes.  A  plane  can  be  represented  by  the  equation 

Ax  +  By  +  Cz  +  D  =  0  or  in  vector  notation  A-r  +  D  =  0 


where  the  vector  A  is  directed  along  the  normal  to  the  plane  with  r  being  the 
position  vector.  If  the  z  axis  is  located  parallel  to  the  wave  wires  and 
perpendicular  to  the  cross  structure  at  the  gage  location  then  the  directional 
cosine  of  the  normal  relative  to  the  z  axis  also  defines  the  angle  of  the  free 
surface  with  respect  to  the  cross-structure.  In  this  case  the  z  coordinate 
represents  the  free  surface  relative  to  the  cross-structure.  The  appropriate 
directional  cosine  is  given  by 

~C  sign(D) 


cos  CLz  = 


2  2  2 

'A  +  B  +  C 
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The  subscripted  z  values  correspond  to  the  three  relative  wave  height 
measurements  at  the  instant  of  time  that  the  impact  angle  is  calculated  and 
the  x, y  values  locate  the  corresponding  probes  on  the  cross-structure.  As  in 
the  case  of  the  relative  velocity,  the  impact  angle  was  computed  at  the 
instance  just  prior  to  the  wetting  of  the  cros3-structure . 


PRESENTATION  OF  RESULTS 


Plots  of  the  impact  pressure  variation  with  relative  velocity  are 
presented  in  Figures  4  and  5  respectively  for  zero  3peed  and  8  knots 
(equivalent  full  scale)  in  head  3eas. 


Figure  4  - 
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Figure  5  -  Impact  Pressure  Variation  with  Relative  Velocity 
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Figure  7  -  Impact  Pressure  Variation  with  Impact  Angle  - 
8  Knots  Ship  Speed 

The  pressure  data  do  not  show  any  consistent  trend  with  variation  in 
impact  angle;  although,  the  zero  speed  data  presented  in  Figure  6,  show 
slightly  less  scatter  than  the  results  obtained  for  8  knots  shown  in  Figure  7. 
A  contributing  factor  causing  the  scatter  may  be  the  inability  to  determine 
the  velocity  with  sufficient  accuracy  at  the  instance  of  impact.  Numerical 
differentiation  is  inherently  difficult  because  of  its  tendency  to  magnify 
small  discrepancies.  Also,  at  the  time  of  impact  the  relative  motion  signal 
becomes  saturated  and  no  longer  indicative  of  the  relative  motion. 
Consequently,  the  relative  impact  angle,  as  well  as  the  relative  velocity,  are 
taken  as  the  values  that  occur  just  slightly  before  impact.  Another  possible 
source  of  discrepancy  may  due  to  the  large  panel  area  over  which  the  pressure 
is  measured  and  its  dynamic  characteristics.  In  addition  to  these  errors 
there  appears  to  be  other  factors,  not  presently  accounted  for,  that  influence 
the  magnitude  of  the  impact  pressure  experienced  by  the  ship. 

Since  there  does  not  appear  to  be  a  clear  relationship  between  the  impact 
pressure  and  the  relative  velocity  and  angle  of  impact  it  was  necessary  to 
examine  a  slightly  more  general  approach.  Dinsenbacher  has  demonstrated  that 
the  impact  pressures  can  be  describe  statistically  by  an  exponential 
probability  distribution  function  where  the  probability  P  of  the  pressure  p 


being  les3  than  some  value  "a"  is  given  by 


c 


where,  E  -  mean  square  pressure. 

and,  C  -  constant  -  ^2 . 

Figure  8  shows  that  equation  [1]  is  fairly  representative  of  the  two  sets 
of  data  presently  being  examined. 


Figure  8  -  Impact  Pressure  Probability 


Ochi  and  Motter^  had  previously  obtained  similar  results  for  bottom  slamming 


by  assuming  that  the  relative  velocity  is  a  narrow  band  Gaussian  process  and 


he  impact  pressure  is  proportional  to  relative  velocity  squared 


where. 


p  -  2kxV2 

is  a  constant. 


The  narrow  band  Gaussian  process  defines  the  probability  distribution  of  the 
velocities  at  impact  as  a  Rayleigh  distribution,  and  in  this  case  the 
statistical  properties  of  the  impact  pressures  can  be  related  to  those  of  the 
relative  velocity,  i.e 


- 2  (p  “  p-> 

f  (P)  - - -  e  [2] 

4Mv 


o 

where,  p  =  impact  pressure  =  2^-V^ 

2 

p*=  threshold  pressure  =  2kj.-v 

av= standard  deviation  of  relative  velocity 
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Figure  9  -  Frequency  Distribution  of  Relative  Velocity  - 


Zero  Ship  Speed 


i  j+h  .*  1 


'I? 

•ft 


*V.'l 


y'll 


u 


VELOCrTY 

Figure  10  -  Frequency  Distribution  of  Relative  Velocity  - 
8  Knots  Ship  Speed 

The  pressure  data  shown  previously  in  Figures  3  and  4  indicates  that  the 
velocity  squared  relationship  does  not  fit  the  data  without  considerable 
scatter.  Also,  the  frequency  distributions  of  the  relative  velocity  at 
impact,  shown  in  Figures  9  and  10  for  the  results  obtained  at  zero  speed  in 
sea  state  8  and  at  8  knots  in  3ea  state  9  respectively,  do  not  appear  to 
closely  match  a  Rayleigh  distribution  as  would  be  expected  for  a  Gaussian 
process.  Some  of  this  discrepancy,  as  indicated  previously,  may  be  attributed 
to  the  inherent  difficulty  in  measuring  the  velocity. 

In  spite  of  the  above  lack  of  strong  correlation  between  velocity  and 
pressure  an  attempt  has  been  made  to  establish  a  statistical  relationship 
between  impact  pressure  and  the  relative  velocity.  If  it  is  assumed  that  the 
root  mean  squared  pressure  is  proportional  to  the  mean  squared  relative 
velocity 


C  V  -co-2 

O  0  V 


where , 


p^  =  mean  squared  pressure 
CQ  -  arbitrary  constant 

Ov»  standard  deviation  of  relative  velocity 

~2 

V  =  mean  squared  relative  velocity 


'Aw 


Y»\*y 
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then,  Dinsenbacher '  s  formulation  [1]  and  Ochi  '  5  formulation  [2]  are 
interchangeable  when  the  threshold  pressure  or  velocity  is  zero.  Reference 
(5)  has  demonstrated  that  the  relative  motions,  hence  Ov,  can  in  most  cases  be 
accurately  estimated  using  the  SWATH  Seakeeping  Assessment  Program  (SSEP) 
avoiding  the  need  for  experimentally  determining  rms  pressure . 

In  order  to  partially  verify  this  approach,  calculations  were  made  of 
some  of  the  basic  relationships  predicted  as  a  consequence  of  the  above 
assumptions  and  compared  with  those  obtained  experimentally.  These  results 
are  presented  in  Table  1.  First,  the  average  number  of  level  crossings  (N) 
were  calculated  for  the  zero  level  and  the  cross-structure  level  at  the  gage 
location  using  the  formulation  for  a  stationary  Gaussian  process 


where,  H  =  displacement  at  level  crossing 

dr  =  standard  deviation  of  displacement 
CTV  =  standard  deviation  of  velocity. 

These  results  were  compared  with  the  actual  number  observed  which  is  shown  in 
Tabie  1. 

Table  1  shows  that  the  calculated  estimates  are  in  very  good  agreement 
with  the  observed  results.  It  was  also  observed  that  the  number  of  impacts 
was  very  close  to  the  number  of  level  crossings  and  that  a  threshold  velocity 
analogous  to  that  found  for  monohull  bottom  slamming  does  not  appear  to  be  a 
significant  factor  in  the  cros3-structure  slamming  of  the  SWATH  ship.  Table  1 
also  presents  the  constant  CQ  computed  by  taking  the  ratio  of  the  root  mean 
square  pressure  to  the  standard  deviation  of  the  relative  velocity  derived 
from  its  spectrum. 


TABLE  1  -  Level  Crossings  and  Other  Impact  Parameters 


8  KNOTS  -  SEA  S~ATE  9 

MEASURED 

COMPUTED  FROM 

SPECTRUM 

MEASURED 

* 

NO.  OF  ZERO  CROSSINGS 

0.46 

0.43 

0.54 

3.52 

NO.  OF  LEVEL  CROSSINGS 

0.13 

0.12 

mm 

0.2" 

RATIO  OF  RMS  PRESSURE  TO 

MEAN  SQUARED  RELATIVE  VELOCITY 
TIMES  DENSITY 

16 

26 

*  Crossings  per  second  with  positive  slope  (model  scale) 

It  may  be  noted  that  the  value  of  this  coefficient  obtained  at  8  knots  in  sea 
state  9  is  slightly  higher  than  the  value  for  zero  speed  in  sea  3tate  8.  The 
two  data  spots  are  insufficient  to  determine  if  the  difference  reflects  a 
genuine  trend  or  simply  experimental  variations. 

If  equations  [3]  and  [4]  are  assumed  to  apply  to  the  SWATH 
cross-structure  slamming  then  according  to  Ochi's  formulation  the  extreme 
value  of  impact  pressure  is  given  by 

Pn  <«)  =  -  lr.  {  i  _  (  i  -  a  )  }  [5] 

h 

where,  t--n  ( ot)  is  the  extreme  value  of  the  impact  pressure  whose  probability  of 
being  exceeded  in  n  impacts  is  equal  to  a.  Equation  [5]  is  only  applicable  to 
operating  conditions  in  which  the  probability  distribution  and  associated  rms 
impact  pressure  are  the  sane,  i.e.  same  sea  state,  heading  angle  and  speed. 
An  equivalent  expression  can  be  found  for  assessing  the  probability  of  the 
extreme  value  in  a  combination  of  discrete  sea,  heading  and  speed  conditions 
provided  each  of  the  probability  distributions  are  known.  Assuming  that  the 
probability  density  of  the  pressure  in  each  condition  is  exponential,  the 
probability  a  of  exceeding  an  extreme  value  p  is  given  by 


a  =  i  -  rr  ( 


l  -  e 


p 


r  c  t 

L  0  J 


where  n- 
and,  "VE, 


number  of  slams  associated  with  i1'*’1  condition 
rms  impact  pressure  associated  with  i*"*"1  condition. 


Equation  [6]  can  be  solved  numerically  to 
associated  with  a  particular  probability 
associated  with  an  assumed  extreme  pressure. 


determine  an  extreme  pressure 
or  conversely  the  probability 


SUMMARY  AND  CONCLUSIONS 

Measurements  of  SWATH  cross-structure  impact  pressures  and  the  relative 
motions  between  the  free  surface  and  the  structure  at  the  point  of  impact  on  a 
1/22  scale  model  of  the  T-AGOS  were  analyzed  to  determine  the  fundamental 
relationship  between  these  parameters  which  could  lead  to  a  procedure  for 
estimating  the  pressure  from  ship  motion  computations.  The  impact  pressure 
was  anticipated  to  vary  as  the  square  of  the  relative  velocity,  but  the  data 
showed  considerable  scatter  when  plotted  in  this  fashion.  Taking  into  account 
the  angle  of  the  free  surface  at  the  point  of  impact  did  not  improve  the 
correlation  as  was  anticipated.  The  above  disparity  can  in  part  be  attributed 
to  the  inability  to  accurately  determine  relative  velocity  which  required 
numerical  differentiation  of  the  displacement  and  possibly  due  to  the  large 
panel  area  over  which  the  pressure  was  averaged. 

Computation  of  the  number  of  level  crossings  at  the  point  of  impact  from 
measurements  of  the  power  spectrum  agreed  closely  to  the  actual  number 
observed.  There  was  little  difference  between  the  number  of  level  crossings 
and  the  number  of  impacts;  therefore,  the  number  of  slams  can  be  assumed  equal 
to  the  number  of  level  crossings  with  little  error.  This  would  suggest  that 
the  number  of  level  crossings  at  the  point  of  impact  can  be  used  directly  as 
a  crude  assessment  of  the  slam  characteristics  of  the  SWATH  3hip. 
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